Smoking-related dysregulation of plasma circulating microRNAs: the Rotterdam study.

BACKGROUND: MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression. Differential miRNA expression, which is widely shown to be associated with the pathogenesis of various diseases, can be influenced by lifestyle factors, including smoking. This study aimed to investigate the plasma miRNA signature of smoking habits, the potential effect of smoking cessation on miRNA levels, and relate the findings with lung cancer incidence. RESULTS: A targeted RNA-sequencing approach measured plasma miRNA levels in 2686 participants from the population-based Rotterdam study cohort. The association between cigarette smoking (current versus never) and 591 well-expressed miRNAs was assessed via adjusted linear regression models, identifying 41 smoking-associated miRNAs that passed the Bonferroni-corrected threshold (P < 0.05/591 = 8.46 × 10-5). Moreover, we found 42 miRNAs with a significant association (P < 8.46 × 10-5) between current (reference group) and former smokers. Then, we used adjusted linear regression models to explore the effect of smoking cessation time on miRNA expression levels. The expression levels of two miRNAs were significantly different within 5 years of cessation (P < 0.05/41 = 1.22 × 10-3) from current smokers, while for cessation time between 5 and 15 years we found 19 miRNAs to be significantly different from current smokers, and finally, 38 miRNAs were significantly different after more than 15 years of cessation time (P < 1.22 × 10-3). These results imply the reversibility of the smoking effect on plasma levels of at least 38 out of the 41 smoking-miRNAs following smoking cessation. Next, we found 8 out of the 41 smoking-related miRNAs to be nominally associated (P < 0.05) with the incidence of lung cancer. CONCLUSIONS: This study demonstrates smoking-related dysregulation of plasma miRNAs, which might have a potential for reversibility when comparing different smoking cessation groups. The identified miRNAs are involved in several cancer-related pathways and include 8 miRNAs associated with lung cancer incidence. Our results may lay the groundwork for further investigation of miRNAs as potential mechanism linking smoking, gene expression and cancer.

Targeted DNA methylation analysis and prediction of smoking habits in blood based on massively parallel sequencing.

Tobacco smoking is a frequent habit sustained by > 1.3 billion people in 2020 and the leading preventable factor for health risk and premature mortality worldwide. In the forensic context, predicting smoking habits from biological samples may allow broadening DNA phenotyping. In this study, we aimed to implement previously published smoking habit classification models based on blood DNA methylation at 13 CpGs. First, we developed a matching lab tool based on bisulfite conversion and multiplex PCR followed by amplification-free library preparation and targeted paired-end massively parallel sequencing (MPS). Analysis of six technical duplicates revealed high reproducibility of methylation measurements (Pearson correlation of 0.983). Artificially methylated standards uncovered marker-specific amplification bias, which we corrected via bi-exponential models. We then applied our MPS tool to 232 blood samples from Europeans of a wide age range, of which 90 were current, 71 former and 71 never smokers. On average, we obtained 189,000 reads/sample and 15,000 reads/CpG, without marker drop-out. Methylation distributions per smoking category roughly corresponded to previous microarray analysis, showcasing large inter-individual variation but with technology-driven bias. Methylation at 11 out of 13 smoking-CpGs correlated with daily cigarettes in current smokers, while solely one was weakly correlated with time since cessation in former smokers. Interestingly, eight smoking-CpGs correlated with age, and one displayed weak but significant sex-associated methylation differences. Using bias-uncorrected MPS data, smoking habits were relatively accurately predicted using both two- (current/non-current) and three- (never/former/current) category model, but bias correction resulted in worse prediction performance for both models. Finally, to account for technology-driven variation, we built new, joint models with inter-technology corrections, which resulted in improved prediction results for both models, with or without PCR bias correction (e.g. MPS cross-validation F1-score > 0.8; 2-categories). Overall, our novel assay takes us one step closer towards the forensic application of viable smoking habit prediction from blood traces. However, future research is needed towards forensically validating the assay, especially in terms of sensitivity. We also need to further shed light on the employed biomarkers, particularly on the mechanistics, tissue specificity and putative confounders of smoking epigenetic signatures.

Plasma MicroRNA Signature of Alcohol Consumption: The Rotterdam Study.

BACKGROUND: MicroRNAs (miRNAs) represent a class of noncoding RNAs that regulate gene expression and are implicated in the pathogenesis of different diseases. Alcohol consumption might affect the expression of miRNAs, which in turn could play a role in risk of diseases. OBJECTIVES: We investigated whether plasma concentrations of miRNAs are altered by alcohol consumption. Given the existing evidence showing the link between alcohol and liver diseases, we further explored the extent to which these associations are mediated by miRNAs. METHODS: Profiling of plasma miRNAs was conducted using the HTG EdgeSeq miRNA Whole Transcriptome Assay in 1933 participants of the Rotterdam Study. Linear regression was implemented to explore the link between alcohol consumption (glasses/d) and miRNA concentrations, adjusted for age, sex, cohort, BMI, and smoking. Sensitivity analysis for alcohol categories (nondrinkers, light drinkers, and heavy drinkers) was performed, where light drinkers corresponded to 0-2 glasses/d in men and 0-1 glasses/d in women, and heavy drinkers to >2 glasses/d in men and >1 glass/d in women. Moreover, we utilized the alcohol-associated miRNAs to explore their potential mediatory role between alcohol consumption and liver-related traits. Finally, we retrieved putative target genes of identified miRNAs to gain an understanding of the molecular pathways concerning alcohol consumption. RESULTS: Plasma concentrations of miR-193b-3p, miR-122-5p, miR-3937, and miR-4507 were significantly associated with alcohol consumption surpassing the Bonferroni-corrected P < 8.46 × 10-5. The top significant association was observed for miR-193b-3p (ß = 0.087, P = 2.90 × 10-5). Furthermore, a potential mediatory role of miR-3937 and miR-122-5p was observed between alcohol consumption and liver traits. Pathway analysis of putative target genes revealed involvement in biological regulation and cellular processes. CONCLUSIONS: This study indicates that alcohol consumption is associated with plasma concentrations of 4 miRNAs. We outline a potential mediatory role of 2 alcohol-associated miRNAs (miR-3937 and miR-122-5p), laying the groundwork for further exploration of miRNAs as potential mediators between lifestyle factors and disease development.

Validating biomarkers and models for epigenetic inference of alcohol consumption from blood.

BACKGROUND: Information on long-term alcohol consumption is relevant for medical and public health research, disease therapy, and other areas. Recently, DNA methylation-based inference of alcohol consumption from blood was reported with high accuracy, but these results were based on employing the same dataset for model training and testing, which can lead to accuracy overestimation. Moreover, only subsets of alcohol consumption categories were used, which makes it impossible to extrapolate such models to the general population. By using data from eight population-based European cohorts (N = 4677), we internally and externally validated the previously reported biomarkers and models for epigenetic inference of alcohol consumption from blood and developed new models comprising all data from all categories. RESULTS: By employing data from six European cohorts (N = 2883), we empirically tested the reproducibility of the previously suggested biomarkers and prediction models via ten-fold internal cross-validation. In contrast to previous findings, all seven models based on 144-CpGs yielded lower mean AUCs compared to the models with less CpGs. For instance, the 144-CpG heavy versus non-drinkers model gave an AUC of 0.78 ± 0.06, while the 5 and 23 CpG models achieved 0.83 ± 0.05, respectively. The transportability of the models was empirically tested via external validation in three independent European cohorts (N = 1794), revealing high AUC variance between datasets within models. For instance, the 144-CpG heavy versus non-drinkers model yielded AUCs ranging from 0.60 to 0.84 between datasets. The newly developed models that considered data from all categories showed low AUCs but gave low AUC variation in the external validation. For instance, the 144-CpG heavy and at-risk versus light and non-drinkers model achieved AUCs of 0.67 ± 0.02 in the internal cross-validation and 0.61-0.66 in the external validation datasets. CONCLUSIONS: The outcomes of our internal and external validation demonstrate that the previously reported prediction models suffer from both overfitting and accuracy overestimation. Our results show that the previously proposed biomarkers are not yet sufficient for accurate and robust inference of alcohol consumption from blood. Overall, our findings imply that DNA methylation prediction biomarkers and models need to be improved considerably before epigenetic inference of alcohol consumption from blood can be considered for practical applications.

Epigenome-wide association meta-analysis of DNA methylation with coffee and tea consumption.

Coffee and tea are extensively consumed beverages worldwide which have received considerable attention regarding health. Intake of these beverages is consistently linked to, among others, reduced risk of diabetes and liver diseases; however, the mechanisms of action remain elusive. Epigenetics is suggested as a mechanism mediating the effects of dietary and lifestyle factors on disease onset. Here we report the results from epigenome-wide association studies (EWAS) on coffee and tea consumption in 15,789 participants of European and African-American ancestries from 15 cohorts. EWAS meta-analysis of coffee consumption reveals 11 CpGs surpassing the epigenome-wide significance threshold (P-value <1.1×10-7), which annotated to the AHRR, F2RL3, FLJ43663, HDAC4, GFI1 and PHGDH genes. Among them, cg14476101 is significantly associated with expression of the PHGDH and risk of fatty liver disease. Knockdown of PHGDH expression in liver cells shows a correlation with expression levels of genes associated with circulating lipids, suggesting a role of PHGDH in hepatic-lipid metabolism. EWAS meta-analysis on tea consumption reveals no significant association, only two CpGs annotated to CACNA1A and PRDM16 genes show suggestive association (P-value <5.0×10-6). These findings indicate that coffee-associated changes in DNA methylation levels may explain the mechanism of action of coffee consumption in conferring risk of diseases.

Smoking-related changes in DNA methylation and gene expression are associated with cardio-metabolic traits.

BACKGROUND: Tobacco smoking is a well-known modifiable risk factor for many chronic diseases, including cardiovascular disease (CVD). One of the proposed underlying mechanism linking smoking to disease is via epigenetic modifications, which could affect the expression of disease-associated genes. Here, we conducted a three-way association study to identify the relationship between smoking-related changes in DNA methylation and gene expression and their associations with cardio-metabolic traits. RESULTS: We selected 2549 CpG sites and 443 gene expression probes associated with current versus never smokers, from the largest epigenome-wide association study and transcriptome-wide association study to date. We examined three-way associations, including CpG versus gene expression, cardio-metabolic trait versus CpG, and cardio-metabolic trait versus gene expression, in the Rotterdam study. Subsequently, we replicated our findings in The Cooperative Health Research in the Region of Augsburg (KORA) study. After correction for multiple testing, we identified both cis- and trans-expression quantitative trait methylation (eQTM) associations in blood. Specifically, we found 1224 smoking-related CpGs associated with at least one of the 443 gene expression probes, and 200 smoking-related gene expression probes to be associated with at least one of the 2549 CpGs. Out of these, 109 CpGs and 27 genes were associated with at least one cardio-metabolic trait in the Rotterdam Study. We were able to replicate the associations with cardio-metabolic traits of 26 CpGs and 19 genes in the KORA study. Furthermore, we identified a three-way association of triglycerides with two CpGs and two genes (GZMA; CLDND1), and BMI with six CpGs and two genes (PID1; LRRN3). Finally, our results revealed the mediation effect of cg03636183 (F2RL3), cg06096336 (PSMD1), cg13708645 (KDM2B), and cg17287155 (AHRR) within the association between smoking and LRRN3 expression. CONCLUSIONS: Our study indicates that smoking-related changes in DNA methylation and gene expression are associated with cardio-metabolic risk factors. These findings may provide additional insights into the molecular mechanisms linking smoking to the development of CVD.

Multi-Omics Analysis Reveals MicroRNAs Associated With Cardiometabolic Traits.

MicroRNAs (miRNAs) are non-coding RNA molecules that regulate gene expression. Extensive research has explored the role of miRNAs in the risk for type 2 diabetes (T2D) and coronary heart disease (CHD) using single-omics data, but much less by leveraging population-based omics data. Here we aimed to conduct a multi-omics analysis to identify miRNAs associated with cardiometabolic risk factors and diseases. First, we used publicly available summary statistics from large-scale genome-wide association studies to find genetic variants in miRNA-related sequences associated with various cardiometabolic traits, including lipid and obesity-related traits, glycemic indices, blood pressure, and disease prevalence of T2D and CHD. Then, we used DNA methylation and miRNA expression data from participants of the Rotterdam Study to further investigate the link between associated miRNAs and cardiometabolic traits. After correcting for multiple testing, 180 genetic variants annotated to 67 independent miRNAs were associated with the studied traits. Alterations in DNA methylation levels of CpG sites annotated to 38 of these miRNAs were associated with the same trait(s). Moreover, we found that plasma expression levels of 8 of the 67 identified miRNAs were also associated with the same trait. Integrating the results of different omics data showed miR-10b-5p, miR-148a-3p, miR-125b-5p, and miR-100-5p to be strongly linked to lipid traits. Collectively, our multi-omics analysis revealed multiple miRNAs that could be considered as potential biomarkers for early diagnosis and progression of cardiometabolic diseases.

Validated inference of smoking habits from blood with a finite DNA methylation marker set.

Inferring a person's smoking habit and history from blood is relevant for complementing or replacing self-reports in epidemiological and public health research, and for forensic applications. However, a finite DNA methylation marker set and a validated statistical model based on a large dataset are not yet available. Employing 14 epigenome-wide association studies for marker discovery, and using data from six population-based cohorts (N = 3764) for model building, we identified 13 CpGs most suitable for inferring smoking versus non-smoking status from blood with a cumulative Area Under the Curve (AUC) of 0.901. Internal fivefold cross-validation yielded an average AUC of 0.897 ± 0.137, while external model validation in an independent population-based cohort (N = 1608) achieved an AUC of 0.911. These 13 CpGs also provided accurate inference of current (average AUCcrossvalidation 0.925 ± 0.021, AUCexternalvalidation0.914), former (0.766 ± 0.023, 0.699) and never smoking (0.830 ± 0.019, 0.781) status, allowed inferring pack-years in current smokers (10 pack-years 0.800 ± 0.068, 0.796; 15 pack-years 0.767 ± 0.102, 0.752) and inferring smoking cessation time in former smokers (5 years 0.774 ± 0.024, 0.760; 10 years 0.766 ± 0.033, 0.764; 15 years 0.767 ± 0.020, 0.754). Model application to children revealed highly accurate inference of the true non-smoking status (6 years of age: accuracy 0.994, N = 355; 10 years: 0.994, N = 309), suggesting prenatal and passive smoking exposure having no impact on model applications in adults. The finite set of DNA methylation markers allow accurate inference of smoking habit, with comparable accuracy as plasma cotinine use, and smoking history from blood, which we envision becoming useful in epidemiology and public health research, and in medical and forensic applications.